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Abstract 

We study the effects of a non thermal neutralino production, due to the late decay in the early 
universe of a single modulus field, in the context of the deflected anomaly mediated scenario. In the 
regime in which the average number of neutralino produced in each modulus decaying process is 
-^LSP 1 also models with a thermal relic density below WMAP data became acceptable models. 
We find out that these models belong to three different classes with the common feature that the 
low thermal relic density is entirely due to coannihilation effects. The neutralino annihilation cross 
section for these classes of models is not particularly high compared with the highest cross sections 
attainable in the generic framework of the MSSM. Hence the detection prospects either by direct 
or indirect WIMP search experiments are not encouraging. 
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I. INTRODUCTION 



The latest WMAP data [l| determine in a very precise way the matter content of the 
Universe. For the cold dark matter (CDM) component the WMAP result is 

ficDAf/i^ = 0.110 ±0.007 (1) 

where ^cdm = Pcdm/Pc is the ratio between the CDM density Pcdm and the critical 
density of the Universe Pc = 1.8791/i^ ■ 10~^^gcm~^ with h the Hubble constant in units 
of lOOkms"^ pc~^. Eq. ([T]) is able to put a strong constraint on the parameter space of a 
generic supersymmetric model. In fact the lightest supersymmetric particle (LSP), usually 
the neutralino, is one of the best motivated cold dark matter candidate (as long as the R- 
parity is conserved) and one can compute the relic density of the lightest neutralino solving 
the following Boltzmann equation 

^ + 3Hn^ = - {av) (n^ - n%) (2) 

where is the numerical neutralino density, H is the Hubble constant, {av) is the averaged 
thermal cross section times the neutralino pair relative velocity while Ueq is the Boltzmann 
equilibrium number density. The physical picture described by this equation is simple and 
the numerical solution is even straightforward until one takes into account coannihilations 



with partic 



can show [2I that the following relation holds 

Pc 

where is the neutralino mass. With f2thcrmai^2 denote the solution of the standard 
Boltzmann equation ([2]). Typically when one compute eq. ([2]) in the mSUGRA framework, 
in which the LSP is a very pure bino in almost all the parameter space, one finds that the 
relic density easily exceeds the WMAP result j^. The same situation happens in other 
supersymmetry breaking scenarios. Moreover the relation ((31) roughly implies the existence 
of some upper bound on the annihilation cross section. Hence if we want to study models 
with cross section higher than that allowed by the thermal relic density constraint we have 
to consider scenarios in which the condition (JS]) can be relaxed. 

The relic density computation can be heavily affected in a non standard cosmological 
scenario which introduces a non thermal source of cold dark matter production. Among 



es degenerate in mass with the lightest neutralino. In a first approximation one 



the different scenarios tfie neutralino production from late decays of scalar or moduli fields 
is an elegant and "natural" possibility [4, 15|. In particular in supergravity and superstring 
theories, where in general the scalar potential has many flat directions, the presence of 
moduli fields that acquire masses from supersymmetry breaking is a common feature (for 
a review in the string theory context see [6!|). Moduli fields have important cosmological 
consequences although their couplings are suppressed by inverse powers of the Planck mass. 
They dominate the energy density in the early Universe and if the modulus mass is suffi- 
ciently low (for example of the order of the electroweak scale) it can destroy the standard 
big-bang nucleosynthesis scenario. This problem is usually termed as cosmological moduli 



problem (see for example 
this problem was given in 



7| and references therein). One interesting possibility to avoid 
^. In that scenario the modulus mass is setted by the gravitino 



mass which is quite high in anomaly mediated models. In this paper we will study the same 

n n □ 

solution applied in the context of the deflected anomaly mediation [8|, 19|, llOj in which the 
tachyon problem of the minimal anomaly mediation is elegantly solved. In particular we 



concentrate on |10j| where we performed a detailed analysis of the neutralino sector of this 



scenario. 

The outline of the paper is as follows: in sec. [Til we describe how the presence of a 
modulus field changes the standard relic density computation in the framework of a generic 
supersymmetry breaking scenario. In sec. IIIII we apply the results of the previous section to 
the deflected anomaly case. In sec. HV] we study the detection prospects for this scenario. 
The last section is devoted to the conclusions. 



II. RELIC DENSITY AND MODULI FIELDS 



4| with a single modulus field 



We closely follow the analysis of Moroi and Randall 
We consider a system of coupled Boltzmann equations describing the time evolution of the 
wimp number density n^, the modulus number density and the radiation density Prad' 

^ + SHn^ = NLSp'^4,n<t> - {(xv) (n^ - n^J 
+ 3Hn^ = -T^n^ 

^ + 411 Prad = ijn^ - iVLSpm^)r^ra^ + 2m^ {av) {n\ - nig) (4) 
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where is the decay rate of the modulus field and TVlsp is the average number of wimp 
produced in each modulus decaying process. These are the two crucial parameters in this 
description. The modulus field decay rate can be formally written as a sum 



Er? (5) 



where each single decay rate F^*^ is associated to one of the possible interactions terms of 
the lagrangian in which is coupled to one of the fields in the neutralino sector, namely the 
two gauginos B, W and the two higgsinos Hi and H2. The index i runs over all possible 
contributions. In the positive deflected anomaly scenario all the terms described in appendix 
A of ^1 are present and so we decide to use the same useful parametrization 
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(6) 



where is the modulus field mass and Mp ~ 2.4 ■ 10^^ GeV is the reduced 4-dimensional 
Planck mass. The naive expectation for is 

^ m3/2 = (7) 

where 7713/2 is the gravitino mass and is the auxiliary supergravity field vev that breaks 
supersymmetry in the anomaly sector. The last equality holds both in the minimal and 
in the deflected anomaly scenario. This relation guarantees that the modulus mass is in 
general quite high. The presence of the modulus field changes the thermal history of the 
early universe [uj. In particular in order to preserve a successful description of the big-bang 
nucleosynthesis we must require that the reheating temperature associated to the modulus 
field decay 

Tnn - 10 MeV ( ) . (^V . (g) 

^" VlOOTeV/ \Mp ) VlO.75/ ^' 

be of order Trh — 10 MeV. In the previous relation Ae// is the effective interaction scale of 
the modulus and counts the total number of effectively massless degrees of freedom. If 
we impose this condition then ~ 100 TeV, assuming that Ae// — Mp and (7* = 10.75. 
Taking into account ([7]) this implies the following condition on the auxiliary supergravity 
field vev 

> 100 TeV (9) 
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The other condition that has to be satisfied to have a source of non thermal neutralino 
production is 

Trh < Tdec ^ (10) 

where with T^ec we have denoted the neutralino decoupling temperature. By definition this is 
the temperature at which the annihilation rate becomes smaller then the Hubble expansion 
rate. The condition (fTOl) is necessary because in order to have a neutralino production by 
the decay of this must happens far away from the thermal and chemical equilibrium. To 
compute the relic density one must numerically solve the system of equations (l4|) . However 
it is possible to give a simple estimate |4| of the solution in two different regimes 

2(27rV45)9.T^^(a^> Pc/s ^^LSP - 1 

(11) 

3jVLSpm;,r^Mg fe2 - 

where pd s ^ 3.6 • 10~^ ■ h? GeV is the critical density to entropy ratio. A key observation 
about eq. (fTTi l is that in the case of very low number of LSP produced in each decay the 
relic density is completely unrelated to the thermal averaged cross section {crv). 

III. NON THERMAL RELIC DENSITY AND DEFLECTED ANOMALY 

We now apply the results of the previous section to the defiected anomaly mediation case 
and in particular to the scenario proposed in (lo| . This scenario differs from minimal anomaly 
mediation for the presence of an additional gauge mediated supersymmetry breaking sector 
that avoids tachyonic masses for the sleptons (which do not transform under the SU{?i) 
gauge group). In this sector there is a gauge singlet chiral superfield X = {Ax,"^ x, Fx) 
which is directly coupled to Nf copies of messenger chiral superfields $i and $i, transforming 
under the fundamentals and anti-fundamentals of the standard MSSM gauge groups (for a 



comprehensive review see [l2fl). The soft breaking sector of the model is specified, as we 
saw in the previous section, in terms of the auxiliary supergravity field vev that breaks 
supersymmetry in the anomaly sector inducing as well a breaking in the additional gauge 
sector and in terms of the ratio //m, where / = {Fx) is the auxiliary gauge singlet field vev 
and m is the messenger mass scale. 

In [l^ we computed the standard thermal relic density without any non thermal contri- 



bution using DarkSUSY 
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taking into account all the coannihilations. In that case the 



auxiliary supergravity field can be as low as ^ 10^ GeV. The WMAP allowed models 
are those in which the neutralino is a very pure bino or a very pure higgsino depending on 
the number of messenger fields Nf (for Nf > 2 a higgsino branch open up and this branch 
becomes greater for Nf ^ 5). We declare that a neutralino is gaugino-like (in particular 
in our case bino-like) if Zg > 0.9 while is higgsino-like when Zg < 0.1, where Zg is the 
usually defined gaugino fraction. In all the intermediate cases we denote the neutralino as 
mixed-like. 

In the new scenario with a non thermal contribution in order to satisfy the constraint on 
the reheating temperature (see eq. ([8])) we must impose 

> 10^ GeV (12) 

that cuts out all the lower part of the parameter space explored in jl^. The condition f fTOll 
on the decoupling temperature turns out to be always satisfied. Moreover we want to add 
another constraint to select models with the highest possible annihilation cross sections. 
This condition, having in mind the relation ((31), can be translated as an upper bound on the 
thermal relic density 

where 0.09 is the 3cr lower bound of the WMAP result This allows us to use the results 
of the detailed scan performed in [lo| . 

Applying both constraints (fT2l) and (fTSll it is possible to rule out the scenario with a high 
messenger mass scale m > 10^ GeV. In fact in this case there is no higgsino branch and the 
neutralino is a very pure bino with a very low annihilation cross section. For a lower messen- 
ger mass scale m ~ 10^ GeV there are essentially three class of models surviving denoted by 
A, B and C. For m < 10^ GeV almost all the parameter space is excluded either due to the 
presence of tachyons or because there is no electroweak symmetry breaking or because the 
neutralino is not the LSP. The result for m ~ 10^ GeV is conveniently summarized in fig. [T] 
where we show in the plane {fl^^'^^^^^h'^^ m^) the three class of allowed models. The plot was 
achieved applying the two bounds (fT2l ) and (fT3l ) on the entire parameter space. In the first 
class of models (class A) the neutralino is a very pure bino with 200 GeV ^ ^ 400 GeV. 
These models are present for Nf < 5. In the second class of models (class B) the neutralino 
is a very pure higgsino with 400 GeV ^ < 900 GeV. These models are always present as 
long as the higgsino branch is open, i.e. Nf > 2. Class B is the one that resembles more the 
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Figure 1: Allowed models with O^^^^^'/i^ < 0.09 as a function of the neutralino mass mv 

scenario considered in where the neutralino is essentially a wino. The wino annihilates in 
a very efficient way into W^W~ and Z^Z^ (when this channels are kinematically opened) so 
the thermal relic density is usually very low. In an analogous way the higgsino has a quite 
efficient annihilation into W^W~ and although the higgsino cross section is lower 

than the wino one of about a factor 10~^. For class B models the non thermal modulus 
contribution could be extremely important to enhance the relic density. In the last class of 
models (class C) the neutralino has a gaugino fraction Zg > 0.1 so it can be either mixed 
or bino like. In both cases it has a very high mass of order > 1.5 TeV. In table [T] we 
show the mass spectrum for three sample models (one of each kind A, B or C). The first five 
rows define the model. The next three rows contain the masses of the lightest neutralino 
of the next to lightest neutralino m^^ and of the lightest chargino . The next six 
rows show the mass spectrum of the squark and slepton sector: the gluino mass nig, the 
stop mass m^-, an average mass of all the other squarks rrig (they are near degenerate), the 
masses of the left and right selectrons, mg^ and mg^ respectively, and finally the mass of 
the lightest stau m^^. The last two rows show the mass spectrum of the Higgs sector: the 
mass of the lightest Higgs m/jO and that of the CP-odd Higgs m^o (quite degenerate with 
the remaining mass eigenstates). 
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class A 


class B 


class C 


(GeV) 


10^ 


1.2 • 10^ 


2 • 10^ 


Fx/m (GeV) 


1.9 • 10^ 


4.9 • 10^ 


9.3 • 10^ 


Nf 


5 


10 


4 


tan P 


50 


10 


50 


sign{n) 


+1 


+1 


+1 




290 GeV 


705 GeV 


2.1 TeV 




912 GeV 


707 GeV 


2.2 TeV 




912 GeV 


706 GeV 


2.2 TeV 


nig 


4.9 TeV 


23 TeV 


20 TeV 


TJlg 


4.2 TeV 


15 TeV 


17 TeV 


mi 


4 TeV 


14 TeV 


16 TeV 




1.2 TeV 


5.4 TeV 


6 TeV 


men 


430 GeV 


2.5 TeV 


2.9 TeV 




294 GeV 


2.5 TeV 


2.5 TeV 




122 GeV 


123 GeV 


124 GeV 




1.4 TeV 


5.2 TeV 


4.4 TeV 



Table I: Mass spectrum for three benchmark models. See the text for the notation. 



As can be seen the scalar sector tends to be heavy with the exception of sleptons for class 
A models. In this case it is worth noting the degeneracy between the lightest neutralino and 
the the lightest stau f. For models B and C the key feature is the near degeneracy between 
the lightest neutralino and the lightest chargino Wi and the next to lightest neutralino X2- 
The thermal relic density of all these models is shown in fig.[Il As can be seen there is a lower 
bound f2^hermai^2 ^ 0.02. Heucc there are no models with an extremely low annihilation cross 
section and this implies (see eq. (fTTi)) that a relic density in the WMAP range could only be 
obtained in the regime iVLSP ^ 1. In fact, as we will see in more detail in the next section, 
for all models a thermal relic density below WMAP data is entirely due to coannihilation 
effects. In fig. [2] we show the allowed values of iVLSP for each value of the LSP mass m^. As 
quoted in [4] the exact value of iVLSP is extremely model dependent. It relies on the precise 
structure of the operators through which the modulus field can decay. In general the range 
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0.001 
0.0008 
0.0006 

Nlsp 

0.0004 
0.0002 



500 1000 1500 2000 2500 
(GeV) 

Figure 2: Contour plot of Q^h'^ in the plane {m^, Ni^sp)- The red (dark shaded) region corresponds 
to the 3cr WMAP result 0.09 < Q^h'^ < 0.13 while the green (light shaded) region corresponds to 
the 5a WMAP result 0.07 < n^h"^ < 0.14. 

of the allowed values is 

10-=^ < iVLSP < 10-^ (14) 
Such range guarantees that all models fulfill the WMAP constraints at least at a 5a level. 

IV. DETECTION PROSPECTS 

In this section we study the detection prospects for the proposed scenario by means of 
direct and indirect WIMP dark matter experiments. 

As can be seen from table [T] the mass spectrum in this scenario is quite peculiar. For class 
A models the LSP is a very pure bino with a mass in the range 200 GeV ^ < 400 GeV. 
The squark and gluino sectors are around 4 TeV and hence they are essentially decoupled 
from the neutralino sector. Also the sfermions are heavy, in particular the charginos and 
the other neutralinos that are above 1 TeV. The exception is the lightest stau f that is quite 
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degenerate in mass with the lightest neutralino. This implies that the dominant annihilation 
channel is XX ~^ t^t^ through a f exchange and the existence of a very efficient neutralino- 
stau coannihilation channel. This latter property make it possible that the thermal relic 
density is well below the WMAP lower bound. From the detection point of view what 
is important is the neutralino-neutralino total annihilation cross section. In the case in 
hand the dominant annihilation cross section (times velocity v) in r+r^ is only of order 
10~^^ cm'^ s~^. Due to the peculiar shape of the energy spectrum of 7s produced in the 
r+r^ channel, the class A models can be studied by carefully analyzing the 7 ray continuum 
flux produced by 7r° mesons decay coming for example from the galactic center jl^. On 
the other hand the cross sections of the two one-loop processes xX ~^ 77 ^iid xX ~^ Z^l 
are suppressed by a factor 10~^ respect to the r+r~ cross section and so are very hard to 



detect 



16|. The direct detection prospects are as well not very encouraging. For class A 



models the typical neutralino-proton spin independent cross section is 

a^p~10-i=^pb (15) 



well below the current limit of the CDMS experiment [1J| and even below the projected 



sensitivity of the upcoming CDMSII [15|. The indirect detection through cosmic rays is not 
favorite as well. To give an estimate we consider the interstellar positron and antiproton 
fluxes at flxed energies, E = 10 GeV for positrons and E = 1 GeV for antiprotons in or- 
der to avoid uncertainties from solar modulation. Taking into account that all the fluxes 
scale as and assuming for the galactic neutralino distribution a NFW profile and the 

n n 

same propagation model used in p/7| and [18|, respectively for the positrons and the an- 



tiprotons, we obtain that the expected interstellar positron flux at -E = 10 GeV is around 
10~^^ GeV~^ cm^^ sr~^ and the expected interstellar antiproton flux at E = I GeV is 
around 10^^ GeV^^ m^^ s~^ sr~^ which are well below the experimental limits determined, 
for a given background, in Even with the recently launched PAMELA satellite (lol . 

that gives an enhanced resolution and a better spectrum determination compared with the 
balloon experiments, the detection possibilities cannot likely be improved. 

In models belonging to the class B the neutralino is a very pure higgsino with 400 GeV ^ 
^ 900 GeV. From table [T] we see that while all the squarks and sleptons are extremely 
heavy and essentially decoupled from the lightest neutralino sector the lightest chargino and 
the next to lightest neutralino are quite degenerate in mass with the lightest neutralino. 
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The dominant annihilation channels in this case are W^W~ and Z^Z^ with a cross section 
(times velocity v) of about 10^^^ cm'^ s~^. As in the previous case the annihilation cross 
section is that required to obtain a thermal relic density in the WMAP range. The fact that 
^^hormai^2 ^ g gg -g g^^jj-giy fjyg the coanuihilatious with the second lightest neutralino 
and with the lightest chargino. For this class of models the one-loop processes XX ~^ 77 ^ind 
XX ~^ Z^l gives an extremely high cross section, with a value of about 2 ■ 10"^^ cm^ for a 
neutralino mass ~ 400 GeV. This could be within the reach of the GLAST sensitivity {2^. 

From the other side the direct detection is completely out of reach given the neutralino- 
proton spin independent cross section 

a^p ~ 10-" pb (16) 
The cosmic ray signals are expected to be higher respect to the previous class of mod- 



els, but still below the sensitivity of PAMELA (see for a general analysis [21|, 122||). In 
fact for class B models the average interstellar positron flux at = 10 GeV is around 
10^^ GeV^^ cm^^ s^^ sr^^ while the average interstellar antiproton flux at ii^ = 1 GeV is 
around 5 • 10^^ GeV"^ m"^ s^^ sr^^. 

Finally class C models have 1.7 GeV ^ ^2.5 TeV. The neutralino can be bino-like 
or mixed. The dominant annihilation channel is bb for any value of the gaugino fraction Zg. 
The direct detection rate is favorite in the case of a mixed-like neutralino, i.e. Zg = 0.5. 
The neutralino-proton spin independent cross section is 

a^p~10-«pb (17) 



so it could be in the reach of CDMSII [15|| although it could be very difficult to detect such 
a high mass neutralino. Neutralinos with Zg = 0.5 or Zg ~ 0.1 are favorite also in the 
case of cosmic ray detection. The interstellar positron flux at = 10 GeV can be as high 
as 4 ■ 10^^ GeV^^ cm^^ s^^ sr^^ while the antiproton flux at E = 1 GeV can be as high as 



10~^ GeV^^ m^^ s^^ sr^^ (for a mass ~ 2.1 TeV). A comparison with the results of 17| 
and [18 1 shows that the maximum attainable positron and antiproton fluxes are respectively 
a factor 10~^ and 10~^ below the highest generic MSSM fluxes. 
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V. CONCLUSIONS 



We studied the effects of the presence of a single modulus field on the relic density 
computation in the the specific framework of the defiected anomaly mediation. The non 
thermal contribution to the relic density allows to resurrect models with Q^^^^^^^h"^ less than 
the WMAP lower bound. These models belong to three different classes with a definite 
neutralino mass range and composition and with a peculiar mass spectrum. The main 
drawback of this scenario is given by the fact that there are no models with a very high 
neutralino pair annihilation cross section. The low value of the thermal relic density is 
entirely due to very efficient coannihilations with particles degenerate in mass with the 
lightest neutralino. As a consequence the detection prospect of this scenario are not quite 
promising. 

The main result of this analysis is that even invoking a non thermal production mechanism 
does not assure in general the existence of models with high detection rate. The consequences 
of a non thermal production has to be carefully analyzed in each specific setup. 
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